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Abstract 

We show that a quartic contact interaction between charged fcrmions can lead to 
Cooper pairing and a superconducting instabihty in the background of a charged asymp- 
toticaUy Anti-de Sitter black hole. For a massless fermion we obtain the zero mode an- 
alytically and compute the dependence of the critical temperature Tc on the charge of 
the fermion. The instability we find occurs at charges above a critical value, where the 
fermion dispersion relation near the Fermi surface is linear. The critical temperature 
goes to zero as the marginal Fermi liquid is approached, together with the density of 
states at the Fermi surface. Besides the charge, the critical temperature is controlled by 
a four point function of a fermionic operator in the dual strongly coupled field theory. 



1 Introduction 



The first successful microscopic description of superconductivity, BCS theory p], describes 
spontaneous symmetry breaking due to a charged fermion bilinear condensate. In the orig- 
inal theory the fermion pairing is driven by an attractive exchange of low energy phonons. 
More generally, the essential feature is a marginally relevant four point interaction between 
excitations about a Fermi surface [2]. Whether this interaction is generated by phonons or 
otherwise is not crucial. The key fact is rather that in BCS-like theories, superconductivity 
emerges from a conventional free Fermi liquid fixed point. An important challenge facing 
condensed matter theory is to characterise the onset of superconductivity from non-Fermi 
liquid states of matter, such as the 'strange metal' phases of high temperature supercon- 
ductors, e.g. [3][1]. 

Recent developments have shown that charged or rotating black holes can carry Fermi 
surfaces [3 [HI [3 El El EHl CI]- We focus on the charged case in 3+1 dimensional asymp- 
totically Anti-de Sitter (AdS) spacetime in what follows; it seems clear that many features 
of the computation will go through in the rotating case also and, with nonzero fermion 
mass, in asymptotically flat space. Near to the black hole, charged fermions experience a 
background chemical potential. If the charge of the fermions is sufficiently large compared 
to their mass, then at low Hawking temperatures they will build up a Fermi surface. The 
essential physics is the same as that for free fermions in flat space with a chemical potential. 
Computations are complicated by the fact that the background spacetime and electrostatic 
potential are nontrivial and vary on scales of order the Compton wavelength of the fermions. 
In general the Dirac equation cannot be solved in closed form on the whole spacetime. 

It is natural to ask whether black hole Fermi surfaces can have BCS instabilities towards 
superconductivity. In this paper we will add a four fermion contact interaction between the 
charged fermions and compute the quadratic term of the one loop effective action for Cooper 
pairs. We will show that under certain circumstances the quadratic action has negative 
modes, indicating a superconducting instability. The effective action is not local in general 
on the curved spacetime background; this complicates finding e.g. the zero temperature gap. 
However, we have been able to obtain an analytic formula for the critical temperature T^. 
Furthermore, for massless fermions we have found the Fermi surface zero mode analytically, 
allowing explicit results without heavy duty numerical work. The critical temperature is 

T,oc/ie-^^^^^'/^^«-, (1.1) 

where ji is the chemical potential provided by the charged black hole at the AdS boundary, 
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Mp is the energy scale of the four fermion interaction, L is the AdS radius, and A'^cfr. is 
the effective density of states at the Fermi surface. From the perspective of the dual field 
theory, the dimensionless quantity MpL determines the magnitude of a four point fermion 
correlator. Figure |4] shows the fermion charge dependence of Aeg., which is given by 



A^cff. ~ kpvF j V^gii^^Hy , (1-2) 

where kp \s the Fermi momentum, vp \s the Fermi velocity, and is the fermion zero mode 



in the black hole spacetime. The precise formula is given in (6.38) below. 

For asymptotically AdS charged back holes the bulk (free) Fermi surface admits a dual 
interpretation, via the applied holographic correspondence (IS [TSj El [El [E] , as a strongly 
interacting (non-) Fermi liquid in 2+1 dimensions O El El E] . It was understood in [8] that 
for fermions with a relatively low charge compared to their mass, the dispersion relation 
of fermion zero modes near the Fermi surface had a non-Fermi liquid form. Furthermore, 
these modes were broad resonances in the spectral density, rather than sharp quasiparticle 
peaks, as is indeed observed in strange metals [17J. The non- Fermi dispersion was shown to 
lead to, for instance, deviations from the venerable Lifshitz-Kosevich formula for quantum 
oscillations [181 [T9j . The black hole BCS instability we present below therefore has the 
potential to dually describe the non-BCS emergence of superconductivity from a strongly 
interacting non-Fermi liquid. Unfortunately, perhaps, we will find that the superconducting 
instability only occurs at larger values of the fermion charge, where the dispersion relation 
is linear (i.e. Fermi liquid like). This can be traced directly to the vanishing of the density 
of states at the Fermi surface in the non- Fermi liquid cases. It may be possible to evade 
this conclusion via alternate bulk pairing mechanisms with long range interactions. 

Recent works have considered instabilities of charged scalar fields in charged black hole 
backgrounds and the corresponding spontaneous symmetry breaking at low temperatures 
|201 [2T1 I22j . If the charge of the boson is sufficiently large compared to its mass [23] it 
will condense, again in strong analogy to the behaviour of charged bosons in flat space 
with a chemical potential. As with the Cooper pairing instability we have just outlined, 
the dual interpretation is of superconductivity emerging from a strongly interacting non- 
Fermi liquid. One difference is that the boson condensation is classical in the black hole 
background whereas for fermions the effect requires a one loop computation, with an ensuing 
nonlocal (bulk) Landau-Ginzburg action. In the Cooper pairing case of interest here, the 
superconducting order parameter is directly related to a fermionic operator in the dual 
field theory. This may be phenomenologically useful and motivates fermion spectroscopy 
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('ARPES') computations in the superconducting state along the hnes of |24j . in which the 
fermion is chosen to couple in a natural way to the bosonic condensate. 



2 Selfinteracting Dirac fermion 

We consider a charged Dirac fermion with quadratic action 

Sm,^, = j (fx^i{i^T^'D^^-m^^) , (2.1) 

where 

D^ = d^ + ^ujabfj"^ - iqA^ , (2.2) 

Tp = ip^r^-, r"'' = rl^T^l, uJab^, is the spin connection, and r^r^ + r^F'^ = 2g'"', with a mostly 
plus metric. We denote bulk spacetime indices by /i, i/, . . . , abstract tangent space indices 
by a,b, - ■ ■ , and specific tangent space indices by underlines as in T-. Eventually we will 
take the background to be a charged black hole in AdS, but for now the metric and gauge 
field are general. 

The BCS mechanism requires an attractive force between like-charge particles. We 
therefore add the simple contact interaction 

^int = / d^x^{^|;,^^^l^){^^^^l;,) , (2.3) 



where Mp is the mass scale of the interaction, = ir-F-r-r-, and t/jc is the charge 
conjugate fermion 

= ci)'^ , c-^vc = -(r'^)^ . (2.4) 



The interaction (2.3), which also appears in color superconductivity [25], is the relativistic 
generalization of s-wave BCS theory: it couples time-reversed, opposite spin states [26l [27] . 
A similar interaction was considered in a closely related context in [21]. However, this choice 
is not unique. Besides choosing a more general contact term, the pairing mechanism could 
arise from exchange of scalar particles, the attractive channel in a nonabelian gauge theory, 
or perhaps graviton exchange. An attractive interaction per se is not sufficient to generate 



superconductivity, but should be in a 'Cooper channel'. The contact interaction (2.3) is 
simpler than an exchange interaction, and can be considered a toy model for these other 
possibilities which may be more natural from the standpoint of string theory on AdS. 
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3 Effective action for the condensate 



Mimicking the standard procedure in BCS theory, we can now perform a Hubbard-Stratanovich 
decouphng to make the action quadratic in spinors. As usual, there is a choice of channels 
to decouple. Given that we are anticipating a superconducting instability of the Fermi 
surface, we choose to decouple in the Cooper channel. Thus we introduce a charged scalar 
A and write the Lagrangian as 

Ant = V^cF^A + ^r^V^cA* - M||A|2 . (3.1) 

Recall that the fermions anticommute. The equation of motion for A sets 

^ = ^^^'^c, A* = ^^,r5v., (3.2) 

and we recover the original action. 

Now consider the Coleman- Weinberg effective action for A, to look for possible instabil- 
ties. Specifically, we will compute the one loop mass term generated for A upon integrating 
out the fermions. The effective action at quadratic order is 

-2 j d^xd'^x'^/^^/g{x')A{x)A*{x')trG'^{x,x')CT^G{x,x')CT\ (3.3) 

Here G{x, x') is the Euclidean Green's function for the Dirac operator in the gauge field and 
spacetime background, G{x,x') = — (^/^(x)V'(x')). G'^{x,x') is the transpose of the Green's 
function in spin indices, i.e. Gjj(x,x') = Gts{x,x'). To derive this expression we used 
C = G^ = —C^, r^-^ = and [C, T^] = 0. See the representation of the gamma matrices 



in equation (4.8) below. Note also that the interaction term in the Euclidean action is 
minus that in the Lorentzian action. We use Lorentzian gamma matrices throughout. 

We now choose coordinates {u,t,x}, with r Euclidean time, and assume the spacetime 
is translationally invariant along {r, x}. In AdS, the radial coordinate is u and the boundary 
directions are {t,x}. Thus we can Fourier transform 

G{x, x')=TY,j ^(^' A:)e-*-"(---')+*(-"-^") , (3.4) 

where the fermionic Matsubara frequencies at temperature T are 

ujn = TTT{2n + 1) . (3.5) 
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We furthermore restrict to configurations in which the condensate A = A(n) only depends 



on the radial direction. The effective action (3.3) becomes 

si^[A] = Ml ^ J du^\Aiu)\^ + dudu'y/giu)g{u')A{u)A*{u')Fiu,u') , (3.6) 
where V2 is the boundary spatial volume and 

F{u, u') = -2T j T^^tr G^{u, u' , iujn, k)Cr^G{u, u' , -icun, -k)Cr^ . (3.7) 

The next step is to relate the Euclidean Green's functions appearing in (3.7) to real time 
Green's functions. This is a little subtle, although the bottom line is that the boundary 
conditions at the black hole horizon mimic the usual effects of finite temperature field theory. 
In particular, as emphasized in |18l I28j . eigenfunctions and eigenvalues, and hence Green's 
functions, are not analytic functions of iojn- This is because regularity at the Euclidean 
'horizon' u = Uj^ typically requires behaviour of the form 

(n-u+)l^"l/(^^^). (3.8) 

The positive and negative thermal frequencies must therefore be analytically continued sep- 
arately. Analytically continuing, by setting ia;„ — )• z, the Euclidean Green's function from 
the upper imaginary frequency axis yields the retarded Green's function G^{z,k), with 
poles in the lower half frequency plane. Analytic continuation from the lower imaginary 
frequency axis gives the advanced Green's function, G'^(z,fc). This relation between Eu- 
clidean, retarded and advanced Green's functions is a general statement that is particularly 
transparent in the black hole context, as we recall in an appendix. 

The sum over Matsubara frequencies can therefore be rewritten as a contour integral 

^5 



T ^ tr G^{u, u', iun, k)CT^G{u, u' , -k)GT^ 

n 

J dztrG'^iu,u',z,k)GT^Giu,u',-z,-k)GTHanh(^^^ , (3.9) 



where the contour G has a segment in the upper half plane and a segment in the lower half 
plane, each going clockwise around the poles of tanh. The analytically continued function 
G has a branch cut on the real z axis. In the upper half plane, schematically, 

G{z)G{-z) = G^{z)G^{-z) , (3.10) 

where G^'^ are the retarded and advanced Green's functions. This product is analytic in 
the upper half plane. In the lower half plane G^ and G^ are exchanged. On the real axis 
the correlators are related by 

G-^iu, u', n, k) = r^-G^{u', u, n, k)^r^- , (3.11) 



5 



where the transpose in G^^ acts on spin indices. This result fohows easily from the definition 



of the various Green's functions, see the appendix. Deforming the contours in (3.9) onto 
the real axis then gives 

/ ^ tanh ^tr r^G^(n', u, n, kyT'-CT''G'^{u, u' , -n, -k)CT'> . 

(3.12) 

Our objective now is to evaluate these integrals. 



4 The charged AdS black hole 

At this point we will specialize to a planar, charged, asymptotically AdS black hole back- 
ground. This is a solution to Einstein-Maxwell theory 

%,.,=/A/=^(^(B+^)-ji,F^). (4.1) 



t2 / 7 2 

ds^ = — -f{u)dt'^ + ij— +dx^ + dy^\ , A = ^{u)dt , (4.2) 



The black hole background is given by 

^ (-f{u)dt^ + 
with 

The horizon is at n = u+ and the conformal boundary is n = 0. The chemical potential fj, 
of the dual CFT is the boundary value of the Maxwell potential 

$ = m(i-^). (4.4) 

We also introduced the ratio of electric and gravitational couplings 



7^ = 4^ . (4.5) 



In terms of the above quantities, the Hawking temperature of the black hole (and temper- 
ature of the dual field theory) is 



47r 47rn-(_ \ 7^ 
The nonzero components of the spin connection are 



6 



Finally, we adopt the following gamma matrix conventions of [SI [23]. These are useful 
for simplifying the Dirac equation once rotational invariance has been used to consider 
momentum in the x direction without loss of generality. Thus 





^ 








V 





-a^ I y cr2 / \ cr2 

(4.8) 

In this representation the charge conjugation matrix is C = T-T-. This representation is 
not quite the same as in [U [23] because our radial coordinate u is the inverse of their radial 
coordinate. 



5 The retarded bulk Green's function and 



In order to compute the effective action of the condensate (3.6) in the black hole background, 



we need the retarded Green's function. It is the unique solution of 

= --^i6^^\x,x'), (5.1) 

subject to certain boundary conditions discussed below. Transforming into momentum 
space except in the radial direction, 

D(n, k)G^(u, u, n, k) = -^i5iu, u) , (5.2) 

where Z)(il, k) is the radial Dirac operator (including the mass term). The Green's function 
equation is solved in the usual manner by multiplying together solutions of the homogeneous 
equation with a discontinuity across the delta function. For a mode with momentum in the 
X direction, 

V' = e-'^^+^'^^VradialC^^) , (5.3) 

the Dirac equation is 

L>(fi, A:)V'radiai(n) = . (5.4) 



This equation is written out explicitly in (6.2) but here we need only some general properties. 
Writing the wavefunction as 

V'radial(^i) = I ' (5-^) 



2 



the two-component spinors 'i/'i and '02 decouple. The gamma matrices (4.8) were chosen as 
in |8ll23] to make this decoupling manifest. 
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The retarded Green's function satisfies ingoing boundary conditions at the horizon. Near 
the AdS boundary (u = 0), the leading solution of the wave equation is 

The two spinors (1,0) and (0, 1) are eigenspinors of T- with opposite eigenvalues, implying 
that Oq, and ba are canonically conjugate (in a radial Hamiltonian slicing). See e.g. |29j . 
Therefore a boundary condition must be imposed on one, with the other allowed to fluctuate. 
When mL > ^ , we must pick the fluctuating piece to be the normalizable mode proportional 
to (0, 1). More generally, for any m we will choose to impose the boundary condition Oq, = 
on the fluctuating modej^ 



With these boundary conditions, the solution of (5.2) is constructed as follows. For each 
= 1,2, take -0™ to be the solution i 
with Cq, = near the boundary. Then 



a = 1,2, take -0™ to be the solution ingoing at the horizon, and ipa'^^ to be the solution 



^l^{u)7Pa'^^{u') U>U' 

(5.8) 

^a'^^{u)ij'^{u') U<U' 



= Gf e , (5.7) 

where 

with ipa = ^V'a'^^- The Wronskian 1^ is a constant related to the conserved charge current: 

At general frequency and momentum {Q, k}, the solution satisfying ingoing boundary con- 
ditions at the horizon will not satisfy the asymptotic boundary condition, so generically 

^bdy 

we have 



V'cf 7^ V'ci We will normalize the solutions so that near the asymptotic boundary u — t- 



^bdy ^ |^ ^ j ^ . . . ^ ^in ^ J_^3/2-mL |^ ^ j ^ | ^ . . . ^ (^^^q) 

We have introduced the quantity Qce{^,k), which is the retarded Green's function of the 
boundary field theory, similar to the original discussion (for bosons) in [32] . The Dirac 



^For mL < |, there is another possible choice ba — discussed in e.g. fS], and for m — there is a 
continuous choice of boundary conditions discussed in e.g. [301 131) . These more general quantization choices 
will not be considered here. 
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equation (5.4) is a real equation for il^a, so the boundary condition implies ■ip^'^^ is real. By 



evaluating the u-independent Wronskian near the asymptotic boundary we obtain 



(5.11) 



We can see immediately that poles of the boundary and bulk Green's functions occur at 
the (in general complex) quasinormal frequencies of the black hole background, at which 
the mode satisfies both the horizon and asymptotic boundary conditions (see e.g. [33]). 



It follows from the gamma matrices (4.8) that the upper and lower spinor projections 
■01,2 are related by A; — )• —k. Henceforth we drop the subscript and rewrite all quantities in 
terms of the first projection, denoting 



i> = i>i , G = Q 



1 ■ 



(5.12) 



We can now rewrite (3.12) using ( |5.8[ ), (5.11), (4.8), (5.5), and rotational invariance, as 

\k\dk 



F(u,u') 



L6 



ianh—G{^,k)*G{-VL,k) x 



2T' 



(5.13) 



^'"(0, uyip'^'i-n, u')ij'^'^y{n, u)^^'^y{-n, u) u<u' 

All wavefunctions in this expression have momentum kx. We have used the identity 

- t^-Gr{u', u, n, kyr^- = Gr{u, u', n, k^ , (5.14) 



or equivalently, using (3.11), GA{u,u',Q,k) = —Gji{u,u' ,0,, k)* . This last statement can 



be seen directly from (5.8) and the fact that ■ip^'^^ is real and V'™ * = see the appendix. 



Note the modulus sign on \k\ and that the integral is over both positive and negative k. 
Positive and negative momenta should be thought of as corresponding to the Gi and G2 
components of the full Green's function, respectively. 



Our first objective is to perform the frequency and momentum integrals in (5.13) (with 
some regulator). In principle all the quantities appearing in this expression could be ob- 
tained numerically. However, we can do better at low temperatures T <^ fi. It was observed 
numerically in O El [7| , and analytically in [8J , that at these low temperatures poles in the 
retarded Green's function can move close to the origin of the complex frequency plane. 
Specifically, the pole location uji,{k) ~ as /c ~ kp- This is of course the signature of a 
Fermi surface in the bulk geometry, as we might expect for charged fermions in a background 
electrostatic potential. Following the experience of BCS theory, which we are essentially 
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replicating in a curved spacetime background, one might expect that the crucial pairing 
physics occurs close to the Fermi surface. 

Close to the Fermi momentum, we will recall in the next section that [8] 

hi 



k^-n/vF + T^^TuiT. 



+ 



(5.15) 



Here k± = k — kp is the perpendicular distance of the momentum from the Fermi sur- 
face, hi and vp are real constants, and is a zero temperature critical exponent to be 
described below. The function T'^'^F^ is the near horizon Green's function, that will also be 
characterised in the following section and in an appendix. At low temperatures, il/T — t- oo, 



(5.16) 



where /i2 is positive and the phase 9 is such that poles of (5.15) are in the lower half complex 
frequency plane. 

We expect the dominant contribution to be due to the singular locus of the boundary 
Green's function (5.15). Thus we restrict consideration to near the Fermi sur fac«g 

\k\dk \kF\dk^ 



27r 



27r 



(5.17) 



and set A; — )• fej? in the remainder of (5.13). We will check the self-consistency of this 



approximation a posteriori. It is simple to perform the integral, leading to 

h\ tanh : 



F{u,u') 



2kF /•°° dn 



n_ 

2T 

■-Q\ 



(5.18) 



^ 2Q/vF + T^''[Mi^)-MT) 

ip^'^yin, u')'fiij^'^y{-n, u')Tp'''{n, u)^'''{-n, u) u>u' 
ip'"^{n, u')^'''{-n, u')^p^'^y{n, u^ip^^^yi-n, u) u<u' 

All wavefunctions in this expression are evaluated at k = kp. In general, as we recall below, 

there can be multiple Fermi surfaces for a given charge [6]. Their contributions will sum in 

the above formula for F(u,u'). 

^As well as near the Fermi surface, another singular region of the Green's function occurs near the 
boundary of the 'log-periodic' region of |8]. Here the denominator of the zero temperature Green's function 
takes the form: g{n, k)'^ ~ 1 + e-(^+*s)^(*)j7-2^^<'=\ where A and B are real constants and A(fc) as 



ku^ —> g7/\/2 (i.e. when u pure imaginary goes to in (6.7 1 below). By explicitly performing the integral 
and then bounding the k integral of ( |5.13[ ) in the dangerous small fl and small \{k) region, one finds that 
this region does not lead to singular low temperature behaviour. The effect of temperature can be estimated 
by replacing tanh ^ by an IR cutoff in the frequency integral at = T. 
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The form of the Green's function (5.15) used in the integral (5.18) is only correct for 



r2 ^ /i. For our computation to be valid we should therefore make sure that the key 
contribution to the integral comes from sufficiently low frequencies. More precisely, we will 
be interested in frequencies in the range T <^ Q <^ fi. This range is the most IR singular, 
because the tanh ^ in the integral can be replaced by 1, and will be seen to capture the 
universal pairing physics. We now evaluate the integral for these frequencies. 

For 17 ^ T one can take the low temperature limit (5.16) of J-i, (p). The denominator 



in (5.18) becomes 2Vl/vf + 2i /i2e"^'' sin0 fi^'^. The constant /i2 is not dimensionless, but 
rather /12 ~ fi^~'^'^. Therefore T <^ Q <^ fi implies, as emphasised in that the first term 
(~ ft) dominates if z/ > ^ while the second (~ fl'^'^) dominates when 1/ < ^. In both cases, 
the range T <^ Q <^ fi implies that the fermion wavefunctions should be simply evaluated 
at = in the extremal T = black hole background]^ Furthermore, given that we are 
then at = and k = kp (at T = 0), the wavefunction is precisely the Fermi surface zero 
mode, 

^0 ^ ^bdy(^ ^ 0, A; = A;f) = i^'^'in = 0,k = kF) . 



(5.19) 



Let us consider the cases v > \ and < ^ in turn. 



1 . 



V > The leading behaviour of (5.18) in an expansion in T ^ /xrg i$ ^ is 



F{u,u') -- 
where 7_e is Euler's constant. 



-Q- log + 7E + log - 



i;°{u)^i;^iu)ij^{u')^i;^iu') , (5.20) 



u < ^: The leading behaviour of (5.18) in an expansion in T ^ /xrg i$ ^ is now 



F{u,u) 



smt 



l-2u\ T ) 



where the constant appearing in the second term is 

Cy = 2(2^'" - l)r(l - 2zy)C(l - 1v 



+ c,)^0(n)VO(n)V'°(n')V°(^^'), 

(5.21) 



(5.22) 



In these expressions /^rg is a renormalisation scale which we implement as a hard cutoff 
on the frequency integral. In any case what will be most important is the temperature 



dependence and, in the former case (5.20), the coefficient of the logarithm. 

^Setting T = r2 = and keeping the coordinates of the wavefunctions it, u' finite means that we miss the 
contribution to the frequency integral from e.g. l — ujuj^ <CT^ One should worry about the possibility 
of IR singular temperature dependence arising from this very near horizon region. In an appendix we check 
that this region does not give additional low temperature divergences to those discussed in the main text. 
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The first of the two cases above is particularly interesting. There is a logarithmic 



divergence associated with low temperatures. The expression in (5.20) is essentially the 
same as that appearing in BCS theory. It is this first case, > ^, in which we will be able 
to consistently describe the onset of superconductivity. The logarithmic divergence in this 
case indicates that the range of frequencies and momenta we have integrated over do indeed 
pick out the dominant contribution to the integrals at low temperatures. 

In the second case there is no low temperature divergence. Recall that the origin of 
superconductivity is the marginal relevancy of the BCS coupling about the Fermi surface [2]. 
What has happenned in the second case is that the modified (non-Fermi liquid) low energy 
dispersion relation, Cl'^'^ ~ k± as opposed to ~ has resulted in the quartic coupling 
becoming irrelevant and hence not leading to strong IR effects. This is familiar from the 



bosonic case [34j. The absence of a low T divergence means that (5.21) is not a controlled 
approximation to the integral and one should consider the full problem before concluding 
beyond doubt that there is no superconductivity in this case. We have shown however that 
there is no instability near the Fermi surface. Consistent with this observation we will find 
below that Tc as = ^ is approached from above. In some contexts, see e.g. 
for a sampling including a color superconductivity case, long range pairing interactions lead 



to additional inverse frequency dependence in the analogue of (5.18). Such interactions 
can compensate for the weaker frequency dependence of the fermion propagator, allowing 
for non-Fermi liquid pairing to occur. In the holographic models under consideration, a 
nontrivial fermion propagator can result from classical propagation on a curved spacetime, 
but the leading fermion interactions are mediated by a contact interaction that seems not 
to allow the non- Fermi liquids to pair. 



In both of the above expressions, (5.20) and (5.21), the u and u' dependence of F{u,u') 
has factorised. This allows us to solve explicitly for the unstable mode and the critical 
temperature Tc. The critical temperature is defined by the appearance of a zero mode for 



the condensate. From the quadratic effective action (3.6) this requires 

M|A(n) + j duy/-g{u')A{u)F{u, u) = . (5.23) 
The factorisation of F{u' , u) immediately allows us to conclude that the critical zero mode 



A°(u) oc /(n)V(^i) , (5.24) 



which seems rather natural. Plugging back into (5.23) leads to a simple formula for the 
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critical temperature (for u > ^) 

M'f = (log ^+7E + log ^) / dnV^gii^'W . (5.25) 



Solving for the critical temperature gives 

Te = f e^-/iRGe-*^^^'/^- (^'^) , (5.26) 
where the effective density of states at the Fermi surface A'gfj (jq) is a dimensionless function 



of the fermion charge q in units of 7, as defined in (4.5). The exponent is the most important 
part of this expression. In the remainder of the paper we will determine this dependence 
of Tc on jq, which is a free parameter in our theory. For our various approximations to 



be reliable we need Tc ^ AiRG- From (5.26) this is seen to hold when ^if-i^ff <■ 1- This 



last inequality can also be thought of as the condition for validity of perturbation theory 



in the quartic fermion interaction of our theory (2.3). That said, the computation of Tc 
is balancing a classical and one loop mass term at a nonperturbatively (in the coupling) 
low temperature. Around and below this temperature a marginally relevant coupling is 
becoming strong and one might worry about the need to resum large logarithms at higher 
orders in perturbation theory. The one loop computation we have performed is in fact the 
only fermion loop contribution to the quadratic effective action for A within our theory and 
so our computation is exact close to the critical temperature. More generally, moving away 
from the quadratic level, in BCS-Eliashberg theory the kinematics of the Fermi surface leads 
to the absence of further relevant operators and the one loop computation is exact, see e.g. 
[2]. It is likely that a similar statement holds for our setup. 

As well as the 75 dependence in the exponent, there is the dependence on j^tj2- The 
Fermi mass Mp does not correspond to a dimensionful scale in the dual field theory (the 
only scales in the otherwise conformal field theory are T and /x). We will discuss in the 
final section how ^.^l is instead related to a dimensionless four point correlator in the 
field theory. This correlator therefore controls Tc in our (dual) strongly interacting field 
theory in the same way the electron-glue coupling controls the critical temperature in a 
perturbative BCS treatment. 



In the prefactor in (5.26) the overall scale is set by ^rg- This is a renormalisation rather 
than physical scale, highlighting the fact that the prefactor is not well defined but is scheme 
dependent. Presumably Tc ~ /i, as /U is the only scale in the theory at low temperatures. 
Renormalisation ambiguities will cancel in dimensionless quantities such as the ratio of 
Tc to the zero temperature mass gap. These ambiguities do not affect the coefficient of 
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the logarithmic divergence at low temperatures, due to frequencies T ^ ^ /i, but will 



shift the order one term in (5.20). Our main interest will be the function NcsXiq) ™ the 
exponent, which is robust. We have primarily kept the prefactor f e"^^ ~ 1.13 in (5.26) in 



order to emphasize the strong similarity with the standard BCS result. The only difference 



between (5.25) and the analogous expression appearing in flat space BCS theory is the 



integral over a spatially dependent zero mode determining the effective density of states at 
the dual field theory Fermi surface. 

As noted above, for large enough q there will be multiple Fermi surfaces. The condensate 



(5.24) and critical temperature (5.26) were given for the case of a single Fermi surface, but 



can be readily generalized to include the contribution from all the surfaces at once. Each 



Fermi surface contributes a term of the form ( 5.20 ) to F(u,u'). Labeling the Fermi surfaces 
with ly > ^ hy n = 1,2,3, ... , with corresponding zero modes '0^^^(n), the solution of the 



integral equation (5.23) has the form 



(5.27) 



Plugging this ansatz into the integral equation gives an eigenvalue equation for Tc. This 
can be written as 



det 



with 



Br, 



Cn 



T\7-total Ty 



Cn 



0. 



(5.28) 



/=5<m)(^)V?„)(^)<)(^)VS'„)(n) . 



The critical temperature is then given by 

vr 



(5.29) 



where N^'^^^ is the largest eigenvalue of (5.28). The corresponding eigenvector is Un, the 



relative contribution of each zero mode to A'^(n) 



We now turn to the computation of the various quantities appearing in (5.25). We have 



been able to find the zero mode analytically when m = 0, allowing many explicit results. 



6 Solution of the massless Dirac equation 

Consider a massive, charged fermion in the black hole background with the wavefunction 



ijj = e' 



-iQt+ikx 



{X+ + X-,i{x+ - X-),0,0) 



(6.1) 
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The combinations x± chosen for convenience. As discussed above, the lower two com- 



ponents of the spinor decouple and so can be set to zero. The Dirac equation (5.4) is 

Lm 



where 



3/ , /' 



3.2) 



D± = fdu 



(6.3) 



The Dirac equation is actually real, and the i appears in (6.2) only because it was inserted 



by hand in (6.1). Generally, this equation must be solved numerically. We will treat the 



massless case, where we have been able to obtain the wavefunctions analytically for small 0. 



Setting 171 = 0, the components of (|6.2|) can be decoupled to give the second order equations 

(6.4) 



where 



V± 



1 

7 



n + qAt± 



if 



4ti2 2u 4 n_i_ 



(6.5) 



6.1 Zero modes at zero temperature 



Our first objective is to find the fermion wavefunctions at = in the T = background. 
From this point on we restrict to the massless m = case as we are able to find an analytic 
solution here. Define dimensionless quantities and radial coordinate z by 



k = ku^ , u = u+(l — z) 



3.6) 



Recall from (4.6) that at zero temperature n+ = ^/Sj/fi. 



Near the horizon, z ~ 0, the wavefunctions behave as z 2 '^'^ with 

1 



We will later be interested in a certain k = kp and define 



3.7) 



(6.8) 



This is the v we referred to in the previous section. We will mainly be interested in cases 
with > 2- We require the z~^^'^ behavior for regularity. 

At zero temperature, the function that appears in the metric can be written as 



/ = 2,z^{z - zo){z - Zq) 



3.9) 
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Plugging this into the decoupled equations (6.4) with $7 = gives the equation for the 



fermion zero modes. It is found to have the exact solution 



V2. 



Z- Zq 
Z- Zq 



X 2-^1 ^ + i^fe ± — g, Z^fc ± 1 + 2z^fc 



l{-l±V2q/zo) 



-2i\/2z 



(6.10) 



3 6 ' ' 3^0(2: — Zq) j 

with N± a normalization. The second independent solution is obtained by replacing Vk 



-Uk in (6.10) 



r]1 = N±x [j^ with i^k 



(6.11) 



with a new normalization A'^-i-. The first solution (6.10) has the required regular behavior 



at the horizon for the 17 = solution. The second will also be required when we consider 



a small nonzero frequency. Inserting the solution into the first order Dirac equation (6.2) 
gives the relative normalizations 



6it^fc - q I Zq 



Zq) 



(6.12) 



Thus we have obtained the zero modes tp appearing in the expression (5.25) for the critical 



temperature. While computing hi and vp will require moving to small frequencies, we 
already have enough information to obtain kp- 



6.2 Fermi momentum kp 

The Fermi surface is characterized by a zero mode that is regular at the horizon and obeys 
certain falloff conditions near the boundary of AdS4. As described in Section [5j the asymp- 
totic boundary condition on the fluctuating mode is 



oc (1 - z) 



3/2 







The equation for the Fermi momentum kp is therefore 

lim(z-l)-3/2(/^ + /_)=0. 



+ • 



2-5>l 



Using (6.10) and (6.12) in (6.14) gives 



Qv + iq 



(6.13) 



(6.14) 



iq 1 
6 ' 2 



v|2,l + 2z.,|(l + iV2 



kp{2i + V2 



(6.15) 
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The hypergeometric functions can be evaluated numerically to solve for the Fermi surface. 
For example, when 75 = 1, we find kp ~ .918528 in agreement with the numerical solution 
of the Dirac equation in [6j. The solutions of (6.15) are plotted in figure [!} For a given q 
there can be multiple Fermi surfaces |6]. When a distinction is necessary, the largest 
will be called the 'first' Fermi surface, the next {kpl the 'second' Fermi surface, and so on. 



10 F 



u+ kp 



-10 







Vk > 






vi<0 - 






Vk > 









qy 

Figure 1: Fermi momentum vs charge of the fermion field. The dashed line is — and the shaded 
region is < J^fc < ^. From left to right are the first, second, etc. Fermi surfaces, which disappear 
when they hit the dashed line. For the vertical axis, recall that m+ = ^/Sj/ii. Positive and negative 
kp correspond to Fermi surfaces in the Green's functions Qi and Q2 respectively. 



We can also use the explicit solution above to make an observation about the Green's 
function at O = 0. The boundary Green's function defined in (5.10), evaluated at zero 
frequency, is 



g{n = 0,k) = i lim 



X- 



1 + x° 



(6.16) 



We will now show that the imaginary part of this Green's function vanishes when is real. 
For imaginary the imaginary part of the Green's function will be positive provided that 
iffc is taken to be positive (this amounts to a choice of sign of a square root in the above). 
For real ffc, as we are interested in, the wavefunctions can be shown to satisfy 

,0 



with P the real phase 

P = TT + 



tan 



X- 



e'^'xl 



2 tan 




\/2g/3) . 



(6.17) 



(6.18) 
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This expression assumes for concreteness that is real and q > 0. It follows that for real 
z/fc the spectral density vanishes away from k = kp, i.e. Im^(a; = 0,k) = 0, while the real 
part has a pole at the Fermi momentum. 



By substituting the Fermi momentum into the zero mode wavefunction (6.10) we can 
plot the radial profile of the unstable Cooper pair mode A*^ of equation (5.24). The result 



for various charges is shown in figure [2] below. The most notable feature of these plots is 
that for large charge the zero modes are supported away from the horizon while as — t- ^, 
at smaller charge, the wavefunctions are supported in the near horizon region. This is 
consistent with the observation, below and in [HI [11] , that for u < ^ the physics of the 
(non-Fermi liquid) Fermi surface is captured by the near horizon AdS2 x region. 




Figure 2: Wavefunctions of the unstable mode A° at different values of the charge for the first Fermi 
surface. From left to right = {11,5,3,2, 1.6}. The curves are normalised so that their maxima 
are equal. At large charge the wavefunctions are supported away from the horizon. As — > | at 
smaller charge {jq — ?> 1.56), the wavefunctions are increasingly supported in the near horizon region. 



6.3 Small frequencies 

Our next goal is to find wavefunctions at T = for small but nonzero frequencies il, using 
a matching procedure. This was done in [8J, although we now have recourse to the exact 



zero mode (6.10). We solve the wave equation analytically in the "near region" defined by 
z <C 1 and in the "far region" defined by z ^ tD (at fixed q). Matching the two solutions 
gives a solution on the full spacetime as long as the near and far regions overlap, 

a}<Cl. (6.19) 



18 



6.3.1 Near region 

In the near horizon region 



f = Gz"^ + 0{z^) . (6.20) 



Similarly expanding all quantities for small z and, anticipating the matching, working to 
lowest order in oj in each power of z, the second order wave equation becomes 

+ 2.xi + ^ {'^f^ T f + 9 - 6P) = , (6.21) 

with solution 

.near _ . , .-^-^..-^ . A , 1 , - . , -, , n.. 



xT' = A±z-2--^e-^^ iFi (^^ ± ^ + i/fc + 1 + 2uk, + B±{uk ^ -Uk) . (6.22) 

On this solution we must also impose the constraint that it satisfies the first order Dirac 
equation. We will do this in the following subsection as part of matching to the asymptotic 
'far' solution. 

Requiring the solution to have ingoing group velocity at the horizon, i.e. leading be- 
haviour of the form e"'"*'^/^^ , fixes 

g^^^^= r(-2..)r(., - f )_ 

and similarly for A^/B^. The ratio = Aj^/B^ is the two-point function in the IR CFT 
living on the boundary of AdS2. In the matching region z ^ Cj, 

^near ^ A±z'l~''>' + 5±z"5+'^* . (6.24) 

6.3.2 Far region 

In the asymptotic 'far' region z ^ uj, we solve the wave equation perturbatively in Co, 
expanding 

xf = X^i^+Cjx^±^ + --- , (6.25) 



X±^ is a combination of the zero modes found in Section (6.1), 



xf=xi+r?i, (6.26) 



where x± rj^ were defined in (6.10), (6.11), and have normalizations A'^-i- and N± 



respectively. Expanding in the matching region z <C 1, 



xf = N^S^z-^-" + N±S^^z-l+'' + 0(w) , (6.27) 
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with 



,^0 



(6.28) 



Comparing the near solution (6.24 ) and far solution ( 6.27 ) in the matching region Co <^ z 



determines the relative contribution of x± and in the far region, 



(6.29) 



Recall that and N- were given in terms of A^4. and 7V+ in (6.12) above. Note that 
obtained in (6.23) is frequency dependent. We will mostly be interested in > ^ and so this 



frequency dependence is subdominant compared to the order uj term in (6.25) and (6.27), 
which we turn to now. 

The first order correction x± satisfies an inhomogeneous wave equation with x± as 
the source. We will need only the leading asymptotic behavior near the boundary z — )• 1. 



The asymptotic behaviour can be found elegantly by integrating the Dirac equation (6.2) 
as in Appendix C of [8] . To leading order at the boundary z — )• 1 one finds 

,3 2z4/;+dnV^(ix°'^ 



1 



\2\ 



L3 



X 



0* 



+ 



(6.30) 



In obtaining this expression there is a contribution at the horizon which vanishes in the 
cases of interest to us (z^ > ^), see [8]. 

6.3.3 Boundary Green's function 

Finally we can obtain the Green's function of the dual field theory living on the boundary 



of AdS4, as defined in equation (5.10). This is given by 



i lim 



X- 



Explicitly in a small Co expansion: 



i lim 



Xa 



^1 X++X- 



X- +??+ -r/° +0(w) 



(6.31) 



(6.32) 



where the wavefunctions x±; V± defined in (6.10, 6.11) are normalized according to (6.12 



6.29), and the asymptotic behavior of the last term in the denominator was given in (6.30). 



Note that due to the matching condition (6.29) 



(6.33) 
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while all other wavefunctions are independent of uj. Near the Fermi surface k = kp -\- k 
the Green's function becomes 

hi 



kj_ -n/vF + h2e^^-^'"'n'^'' ' 



(6.34) 



which is the advertised zero temperature limit of (5.15), using (5.16). In this above formula 

.0 



hi 



VF 



lim 



X- 



L3 / /■"+ 



-1 



lim(l 



zr'ixi 



X 



|2 



(6.35) 
(6.36) 



with all wavefunctions evaluated at k = kp. tp^ is the zero mode spinor with components 



X± as in (6.1). h2,0 are easily obtained from the expression above but will not be needed. 



Note that all of these wavefunctions are hypergeometric functions that have been given an- 
alytically above. From these expressions one finds that hi and vp are real. These quantities 
are plotted as a function of charge in figure |3j Similar plots for the first and second Fermi 
surfaces are found in |8|. 



Vf 0.0 





Figure 3: Left: Fermi velocity vs charge, equation (6.36). vp vanishes at the dashed line when 
ly = Right: hi vs charge, equation (6.35). hi vanishes at the solid line when v = Q. The multiple 



lines in each plot are for the various Fermi surfaces, in ascending order with the first Fermi surface 
on the left. Note that vp and hi have the same sign as kp. As above, positive and negative kp 
correspond to Fermi surfaces in the Green's functions Qi and Q2 respectively. 

The decrease of the Fermi velocity with charge and the fact that the velocity tends to 
the speed of light at large charges in figure [3] suggest a geometrical interpretation. As the 
charge is lowered, the zero mode wavefunction is supported increasingly close to the black 
hole horizon, c.f. figure [2| The gravitational redshift then reduces the local speed of light 
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relative to the boundary value. This observation is made in passing in [8], while similar 
phenomena have been discussed previously in different contexts, e.g. [38] . 



6.4 The critical temperature as a function of charge 



We can now obtain the critical temperature from (5.26), which we rewrite here as 



with 



eff. 



hfvpkp 



vr. 



L4 



(6.37) 



(6.38) 



Using our zero mode (6.10) as well as the above results for {hi,VF, kp} we plot N^s, below 



in figure m The zero mode should be normalised overall according to ( |5.10| ). As well as 
the A'eff. for each Fermi surface individually there is also the total effective density of states 



N^^'^^ and corresponding critical temperature Tc of (5.29). Thus N^g^^ determines the 



actual critical temperature of the system and is also plotted. 

3.0 I 




Figure 4: Lighter lines show the effective density of states at the Fermi surface Ncs. vs charge, 



for each Fermi surface with the first Fermi surface on the left. Recall from (6.371 that Tc 



MpL /Woff. fpjjg dark line is iV^g*^', defined in (5.28), which accounts for the presence of multiple 
Fermi surfaces. 



When the density of states N^s, goes to zero, then clearly from (6.37) the critical tem 



perature Tc also vanishes. This is intuitively reasonable given that it is the states at the 
Fermi surface that pair and condense. Similarly, the larger A'efr. , at fixed MpL'^, the larger 
the critical temperature. Therefore in figure |4] we see that the critical temperature goes to 
zero at jq ~ 1.56, corresponding to = h. Also interesting in this plot is that the fermions 
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that pair at the highest temperature can be on different Fermi surfaces at different charges, 
leading to kinks in Tc as a function of ^q. 

Recah that for our computation to be consistent we need Tc <^ ^. We see in figure [4] 
that for fixed MpL?' this requirement will break down at very large charge. It would be 
interesting to relax this condition by working with general charged black holes rather than 
the near extremal solutions we have considered. In the general case it will likely not be 
possible to solve the equations analytically. 



7 Discussion and dual field theory 

Our main result is figure |4] which shows the critical temperature for the onset of supercon- 
ductivity as a function of the charge of a massless fermion field in a near extremal charged 
black hole background. There are two immediately interesting aspects of this plot. Firstly, 
that at sufficiently large charge there can be crossovers at which different Fermi surfaces 
give the highest critical temperatures. Secondly, that the critical temperature goes to zero 
at the charge such that the critical exponent = with apparently no superconducting 
phase transition from this mechanism for v < ^. 

The regime < 5 contains the non- Fermi liquids described in [8j , including the marginal 
case. The lack of pairing in these cases is potentially disappointing from the point of view 
of dually describing the emergence of nonconventional superconductivity from a non- Fermi 
liquid, and of thereby modeling interesting experimental systems. It is a posteriori not 
surprising however: the vanishing is directly related to a vanishing density of states at 
the Fermi surface which is in turn related to the vanishing of the residue of the Fermi surface 
pole of the Green's function, which is itself the source of the non-Fermi liquid behaviour. 
We note below that one could explore alternative pairing mechanisms to attempt to evade 
this conclusion. Alternatively the onset of superconductivity may not be related to the 
fermions (at least not explicitly), along the lines of |2H I22j. 

Before discussing various loose ends in our above computations and indicating future 
avenues for research, we should characterise the onset of superconductivity from the dual 
field theory perspective. In the bulk we have essentially repeated the BCS computation 
in a curved spacetime. The dual field theory is however strongly coupled. The exponent 



in our formula for the critical temperature (5.26) depends on two parameters, the effective 
density of states at the Fermi surface, N^^^iX-^q)^ and the bulk coupling jjhj^- The effective 
density of states has the same interpretation in the field theory. The bulk coupling however 
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determines the leading contribution to a four-point correlator of the fermionic field. Let 
^ be the (two component) fermionic operator of the dual field theory that is sourced by 
the boundary value of the (two) non-normalisable components of the bulk fermion -0, as 



discussed around equation (5.10) above. Then 

where a particular spin index structure and coordinate dependence follow from the inter- 



action (2.3) and could be computed from Witten diagrams. To summarise: the critical 
temperature is controlled by the magnitude of a specific fermion four point function in the 
zero temperature, zero chemical potential strongly coupled quantum critical theory. This 
fact parallels the observation in [23] that for the bosonic holographic superconductors, the 
critical temperature is determined by the two point function of the bosonic field. In both 
cases the answer is phrased in terms of quantities that are natural at strong coupling. 

In computing the effective action for the Cooper pairs A we have considered the effect 
of virtual fermions, but not of photon or graviton loops. These loops will dress the fermion 
propagator and also renormalise the A-^^ interaction of our decoupled Lagrangian. From 



the Einstein-Maxwell action (4.1) we see that photon and graviton loops will come with 
powers of and j^. Often in string theory realisations 5^ ~ ^ ~ N^"^ ^ 1, where # is 
a positive number and is a measure of the (large) number of degrees of freedom per site 
in the dual field theory. These processes are therefore strongly suppressed compared to the 
term we have computed. There may be interesting physics associated with the resummation 
of these terms which deserves further study. While the long range repulsive photon should 
be screened, the (electric and magnetic) photon 'vertex corrections' may lead to important 
effects if is sufficiently small. To be completely safe with what we have done here we 
could take the strict large limit to set these terms to zero. If we do this we should keep 
j^2^2 finite but small in the large N limit. While this introduces a hierarchy between the 
matter and gravitational sector, such hierarchies are famously known to arise in actually 
existing quantum gravity systems. 

Various future directions present themselves. Perhaps most interestingly, we have not 
discussed a specific pairing mechanism in the black hole background that would generate 



a quartic fermion interaction similar to the one we introduced by hand in equation (2.3). 
While simply introducing a quartic interaction is consistent with a bulk effective field theory 
approach to holography, it may be that various natural candidates for the glue - say scalars, 
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overscreened photons, nonabelian fields or gravitonq^J- lead to distinctive pairing physics 
with interesting dual field theory interpretations. It would be interesting in this context to 
make contact with the physics of colour superconductivity \25\ I39j . Presumably one can 
obtain a condensate with p and d wave pairing. Perhaps different interactions can persuade 
the reluctant modes with u < ^ to condense. To explore pairing of the i' < \ modes, the 
formalism recently developed by [11] may be useful. We noted in our discussion in the main 
text above that in other contexts, e.g. |35 1 136 | [37]. non- Fermi liquids naturally come with 
long range interactions (as opposed to our contact interaction) which can induce pairing 
instabilities in non-Fermi liquids. 

We have also not discussed the physics of the Cooper pair condensate below the critical 
temperature. Natural quantities to compute would be the effective bulk Landau-Ginzburg 
action just below and the zero temperature energy gap. There are two complications in 
doing this relative to the usual flat space BCS gap computation. The first is that when the 
mass of the fermion is of order the curvature of the background spacetime, as is often the 
case in applied holography, then the effective action for the condensate will be nonlocal in 
the radial direction. The second difficulty is that often the most interesting zero temperature 
physics will include the backreaction of the condensate on the geometry. Selfconsistently 
solving simultaneously for the background and the condensate would involve computing 
functional determinants in an unspecified background. This latter complication could be 
avoided by working in a probe limit for the condensate. If the condensate is stabilised by 
higher order terms in the effective action for A (as opposed to via its interaction with the 
Maxwell field, as in [21]) then working in the limit M^L^ ^ {g"^ , -^j should be sufficient. 
In such a probe limit it may be possible to solve the nonlocal gap equation. 

It should be possible to adapt our computation to the case of rotating black holes and 
the Fermi surface discovered in [10]. Furthermore, that paper used an (astrophysically appli- 
cable) WKB approach to reduce the computation of Green's functions to finding geodesies. 
In our context this could appear as e.g. a large mass and charge limit of the bulk fermionic 
fields {v can be kept fixed in this limit). It would be interesting to see to what extent such 
a limit simplifies our one loop calculation and perhaps gives a handle on the question of 
backreaction mentioned above, by giving a local effective action. 
*The possibility of gravitational pairing was mentioned in [24] , 
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A Green's function conventions 



In this appendix we define our conventions for various bulk Green's functions of the Dirac 
fermion ip. The Euclidean, retarded and advanced Green's functions are defined by 

Ge{x,x') = -{TEH^mx')), (A.l) 
GRix,x') = -ie{t-t'){{ij{x),ijix')}), 
Ga{x,x') = ie{t' -t){{^l^{x),^{x')}), 

where T is real time ordering and Te orders with respect to Euclidean time r = it, or more 
generally, — Imt. For real frequency il, the Fourier transformed advanced and retarded 
Green's function are easily seen to be related by 



GAi^, k, u, u) = r^-Gnin, k, u, u)^T^- 



(A.2) 



where the transpose in acts on spin indices. We have used F- = 



t\-i 



The relative sign conventions in ( A.l ) are chosen so that all of the Green's functions satisfy 
the same equation: 

{V^D^ - m)G = i<5(4)/V=5 • (A.3) 
This follows from (F-)^ = — 1 and the equal time commutator 



{Mx),4i^')} = s^'\x,x')/Vh, 



(A.4) 



with h the metric on a spatial slice. 



In the black hole background, using (4.8) and (5.8), acting on one of the two-component 
projections we have 



Therefore from (A.2) 



ia^Gji{uj, k, u , u^ia^ = Gr{oj, k, u, u')* . 



GAi^, k, u, u) = —Gr{uj, k, u, u'y 



(A.5) 



(A.6) 
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In order to confirm that the convention ( A.l ) corresponds to the Ga,r used in the derivation 



of the effective action, note that regularity at the Euchdean horizon in (3.8) uniquely fixes 
the Green's functions that should be used in the Lorentzian formula. For Imz > 0, we must 
pick the ingoing solution, and for Imz < 0, we must pick the outgoing solution. Therefore the 



Green's function G"^ in the derivation of F{u,u') around (3.10) is unambiguously defined 
to be the Green's function in Lorentzian signature with outgoing boundary conditions. 
Therefore Ga is obtained by replacing 



in the retarded Green's function (5.8). From our solution of the Dirac equation, we see that 

Also from the Dirac equation i/j^'^'^ is real, since it has a real boundary condition. Therefore 
we find 

Ga{^, k, u, u') = —Gr{uj, k, u, u'y , (A. 9) 



in agreement with the conventions (A.l) 



B Contributions from the near horizon region 

In the text we have shown that when v > \ there is a logarithmic low temperature divergence 
in the effective mass for the Cooper pairs, as described by the quantity F{u, u') in equations 
(3.6) and (5.13). In evaluating ( |5.13| ) we substituted the zero mode at the Fermi surface 



iIj^{u) for the wavefunctions u). In regions where an expansion in powers of is 

possible, this substitution picks out the most IR singular terms, as higher terms come with 
additional positive powers of 0. However, there are three dimensionful quantities of interest: 
{il, T, n}, and the vanishing limits of these quantities do not commute. In particular, putting 
T = = 0, as we have done in the wavefunction, cuts out the range of integration where 
z < iln+ and also the coordinate range z < Tn+ (recall that u = u+(l— z)). In this appendix 
we check that there are no additional logarithmic (or worse) temperature divergences from 
this near horizon region. To do this we must generalise the discussion in the text to (small) 
finite temperatures. 

In order to study this region, we need the wavefunctions for small Co = Vtuj^ , r = AttTu^ 



and z. The second order radial wave equation in this regime for a spinor of the form (6.1) 
is 

z {6z + t) xl + {I2z + r)x'± + h - ^ ± iq - ^ I6z(6z + r) ) ^± = ° ' ^^'^^ 
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Here q, k and u are as defined in (6.6) and (6.7). The ingoing and outgoing solutions to 
tliis equation are 



T(6z + r)+4 6 ^" T X 



2F1 



1 -r 1 _1_ ,/ 

2 =F 2 + ^fc 



!l 1 -r 1 
6'2^2 



iq 



2iuj 6z 

T T 



(B.2) 



(B.3) 



The normalized 'in' and 'boundary' wavefunctions defined in (5.10) are linear combinations 
of these solutions in the near horizon region 



X± 

bdy 



(B.4) 
(B.5) 



where D,E,F are w, T-dependent normalizations to be determined. 
For the ingoing mode, expanding in the matching region (z ^ ^,T), 



with the finite temperature IR Green's function 



G 



± 



6/ 



;i)r(i 



6 



(B.6) 



(B.7) 



Comparing coefficients to the far region wavefunction (6.27) we find the normalization 



1 



| + ^,)r(i + z., + | 



r(2i/fc)r(i^i 



^)6^2+^^-^i 



(B.8) 



where Sj^ was defined in (6.28), is chosen to normalize the zero mode at the boundary 
such that 



lim(l 



lU, 



3/2 



(B.9) 



and is determined by (6.12). The condition above is the first equality of (5.10). The 
constants A'^-i- characterise the u = T = solution and so do not depend on w and r. This 
fixes the normalization of the 'in' mode. 

For the 'boundary' mode, in the far region we can expand 

V;bdy = ^iii + 0(w,fcx) • (B.IO) 

Therefore the normalizations of ^ to leading order near the Fermi surface, are schemat- 
ically of the form 

F± ~ D±[l + u + k± + {Gi^)-\u + k±)] , 



(B.ll) 
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where all the relative coefficients are independent of Co and have been suppressed. The 
inverse factors of are present so that the Cj dependent term behaves like G:z~^^^ in 
the matching region. Expanded for S> r, the regime of interest shortly, the various 



normalization factors in (|B.4|) and (|B.5|) have frequency dependence 



(B.12) 



B.l Checking for near horizon, low temperature divergences 



Low temperature IR divergences arise in the ^ integral of (5.13) in the range T <C Jl. This 
can be checked explicitly, largely following from the fact that for Vl <^T the tanh term 
gives an extra power of in the numerator. In any case a divergence for < T would 
indicate an IR divergence that is not cured by finite temperature, which is not the physics we 
are after. In the regime T ^ we can replace tanh JJi — )• 1 and T^'^J>(p) — )• in (5.13). 



We will check for IR divergences in the seemingly more dangerous case of > \. Similar 
computations go through for v <\. For the purposes of isolating possible low temperature 



divergences, the integral in (5.13) may therefore be written (with z < z' without loss of 
generality) . 



F{z,z') 



Re 



(B.13) 



We now need to substitute the low temperature, low frequency, near horizon solutions (B.4), 



normalized according to (B.12), into this integral and study the temperature dependence. 



The dangerous regime \s z,z' <t -^Cj because \i t < z, z' <^ uj then the a) integral is cut 
off by z or z' rather than r and there cannot be a temperature dependent IR divergency 



In the regime z, z' < r ^ a) one can expand (B.4) to find e.g. 

"(6z + r)+4 -Q 



(B.14) 



and similarly for x±^ ■ Substituting these expressions into (B.13) leads to various terms 



The fact that < z, z' < r means that powers of z, z' introduce factors of the temperature. 
The only dangerous terms can be seen to take the form 



Fiz.z') 



doj / w z Qz' + T 
cos ( — log — 



(B.15) 



z' uj \2t z' 6z + t 

^One should also worry about 2, — >■ divergences in F{z, 2') that could become temperature dependent 
divergences upon performing the integral in the equation of motion ( |5.23[ ). The most dangerous region is 



T <^ z' < z <^ Cj. Using an argument similar to ( B.16 1 below, we have checked that the small 2, 2' behaviour 
does not lead to new low temperature instabilities. 
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When z = z' this integral scales like logr, while for fixed z ^ z' the zero temperature 
limit is regular due to strong oscillations of the integrand. Thus we find that we do indeed 
have an additional logarithmic low temperature divergence due to the near horizon region. 
However, because this divergence is restricted to the small region < z, z' < r, we can now 
immediately see that it does not affect our results in the main text. In particular, it does 
not lead to a new instability localised in the near horizon region: The equation of motion 



for the zero mode (5.23) requires 

1 



max I Afz) 



max 



dz'F{z,z')/\[z') 



(B.16) 



In the perturbative limit MpL? ^ 1 that we are considering, this equation does not have 
any low temperature solutions. 
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